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Abstract 

Magnetic relaxation in large spin molecular paramagnets is often found 
to behave as 5M(t) ~ \fh at short times t. This behaviour was explained 
by Prokofiev & Stamp as arising from dipole interactions between molecu- 
lar spins. However, as observed by Miyashita & Saito, the same behaviour 
can arise from a different mechanism which, in the present work, is related 
to hyperfine interactions. The Miyashita-Saito scheme is found to be pos- 
sible at short times if the nuclear longitudinal spin relaxation is very slow. 
In the case of moderately slow nuclear spin relaxation, the electronic mag- 
netization variation 5M(t) is initially proportional to t, then to \ft and 
finally to exp(— t/r). This behaviour may be mostly expected in dilute 
systems. 

pacs 05.30.-d, 05.40.Fb, 75.50.Xx 

Synthetic molecular nanomagnetspP provide reproducible microscopic sys- 
tems with a large magnetic moment which may have macroscopic properties. 
The most widely studied materials are, in the usual terminology, Mn^ac (with 
a relaxation time of 2 months for the magnetization at 2 K) and Fes, whose 
faster relaxation allows easier experiments. 

Magnetic relaxation of these materials at low temperature is a challenging 
problem. The material is initially magnetized by a strong magnetic field H~ 
in the easy magnetization direction z. At t = 0, the external field is suddenly 
given the value H ex t, also in the z direction. One measures the magnetization 
M z (t) = M z (0) + 5M z (t). At low enough temperature, the following behaviour 
is observed for short times in Fes [HIS] and Mni2 01 

SM z (t) = 5M z (oo)AVi (1) 

where A is a positive constant. This square root behaviour is in contrast 
with usual relaxation which is exponential, and therefore linear for short times, 
6M z {t) = 5M z {0)A't. 

The square root behaviour was predicted theoretically by Prokofiev & 
Stamp 6^ for the demagnetization of a saturated sample. The molecule of 



interest can be modelled by a 'molecular' spin S, of electronic nature, whose 
modulus s is large (s = 10 for Mni 2 and Feg). This spin is subject to an 
anisotropy hamiltonian which, in Mni2, may be written as a first approximation 

as 

H a = -KS 2 (2) 

where the constant K is positive. 

Relaxation is slow because positive and negative values of S z form two dif- 
ferent potential wells separated by a barrier of height Ks 2 . At temperatures 
and fields of interest (ksT <C Ks 2 and g^BsH ext <C Ks 2 where fis is the 
Bohr magneton and g the Lande factor), thermal activation by phonon ab- 
sorption is not possible and magnetic relaxation takes place by spin tunneling 
through the barrier under the effect of additional terms of the Hamiltonian 
which do not commute with S z . Such a hamiltonian, adequate for Fes, is 
H = Ho + Hi = —KS 2 + giiBH z S z + BS 2 , where the magnetic field H has 
been introduced. 

If H were constant, spin tunneling would be a periodic oscillation of S z 
between —s and s, which could be analyzed by diagonalization of the hamilto- 
nian. In the relaxation mechanism, the environment plays an essential part. It 
will be mimicked for each molecular spin by a time-dependent magnetic field 
H z (t) = H ext + AH Z (t) . The additional component AH Z (t) is produced partly 
by the dipole interaction with the other molecular spins, and partly by the 'hy- 
perfine' interaction with the nuclear spins. The transverse components AH X 
and AHy, which slightly modulate the tunnel splitting 2Tlllit, will be neglected. 
The z-component AH Z , though not larger, is very important. Indeed, spin tun- 
neling is only possible between two eigenstates of J5J which have nearly the 
same energy. This occurs if, and only if the local longitudinal field H z is close 
to or to particular values (or 'resonances') which depend on K . Attention will 
be focussed on one of these values which will be called Hi . Roughly speaking, 
the condition for tunneling is that H z satisfies gnB s \H z — Hx\ < Tilut- 

Prokofiev & Stamp 6 have shown that the dipole interaction with the other 
molecular spins leads to square root relaxation, formula . In their theory, this 
is related to the r -3 behaviour of dipole interactions. It was later suggested by 
Miyashita & Saito[7] that square root relaxation may follow from another, com- 
pletely different mechanism which may be equivalently termed 'Wiener process', 
'random walk', or 'diffusion'. 

In the present note, interactions between molecular spins are ignored. The 
square root law is shown to arise from hyperfine interactions only if nuclear 
spin relaxation is very slow. If it is not so slow, our results are different from 
those of Miyashita & Saito. 

A simplified model will be used. The assumptions are the following. 

1) Interactions between molecular spins are neglected, AH z (t) is the hyper- 
fine field. 

2) The hyperfine field AH z (t) is a sum of independent random components 
H?(t) which can take (2Ii + 1) different values with a definite probability per 
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unit time of jumping from one value to the other. The index i labels the various 
nuclear spins. 

3) The number of nuclear spins interacting with a given molecular spin S is 
large. 

4) Each molecular spin has a probability A per unit time to relax when the 
local field H z {t) = H ext + AH z (t) acting on this spin is comprised between 
Hi — e and Hi + e, where e ~ Twjt I '(d^Bs), otherwise there is no relaxation. 
The value of A is expected to depend on the tunnel splitting, hyperfine field and 
nuclear spin dynamics. 

An extreme case is when a molecular spin has completely relaxed at time t 
if, and only if, the local field H z (t) = H ext + AH z (t) acting on this spin has 
been equal to Hi at some time t' between and t. This limiting case will be 
called 'very slow nuclear relaxation'. 

Assumption (1) is not claimed to be a good approximation, but a simpli- 
fication consistent with our purpose, to investigate y/t relaxation arising from 
hyperfine interactions only. 

Assumption (2) implies that quantum coherence is lost after a short time. It 
is an oversimplification with respect to the real mechanism of nuclear relaxation. 
In Mn 12 below about O.IK, this mechanism seems to be complex, involving 
inhomogeneities of the crystal [S]- 

Assumption (3) is physically realistic since dipole interactions are long ranged. 

Assumption (4) is appropriate when Hi ^ 0, so a spin can tunnel from the 
lowest state of its initial well to an excited state of the other well, where it deex- 
cites with phonon emission. In the case Hi = 0, tunnelling takes place between 
the lowest states of each well and phonons have no effect. Then assumption (4) 
must be reformulated. This will be done in a separate section 



1 Analogy with a random walk 

According to assumption (4) the magnetization at time t depends on the proba- 
bility p(h, t) that the hyperfine field H z (t') has taken the value Hi for < t' < t 
if the initial field was H z (0) = h. The field H z is the sum of contributions H z 
of many nuclei. According to assumption (2), these nuclei flip by random, un- 
corrected jumps. Thus, they are similar to the steps of a random walker In the 
simplest case, the random walker has the same probability to go forward or back- 
ward. Then, the probability 1 p(h,h',t) that a random walker is at AH Z = h' 
at time t if he started from h at time satisfies the diffusion equation 

-p(h,h',t)=D(h')-^p(h,h',t) (3) 

1 To simplify the language, the word 'probability' will often be used instead of 'density of 
probability'. The actual meaning is clear from the formulae. 
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where D is related to the relaxation time of nuclear spins. It will generally be 
assumed to be independent of h! . Then the solution of J2J is 



1 -(h-ti) 2 
p(h,h',t) = —== cxp l " ' (4) 

Formulae © and (0J are not valid for long times. Indeed the field distribu- 
tion g(AH z ) has a finite width AH (Fig. QJ, so that the random walker cannot 
reach fields higher than AH. This can be accounted for by a force f{h'). The 
diffusion equation is thus replaced by the Fokker-Planck equation 



jL p (h,ti,t)-f(ti)p(h,h',t) 



(5) 



which describes an Ornstein-Uhlenbeck process 0. The force is easy to calculate 
by writing the equilibrium condition p(h, h', oo) = g{h'). It follows 

f(h) = ±lng(h) (6) 

If g(h) is a gaussian, g(h) = (2t:)-^ 2 AH- 1 cxp[-h 2 / (2AH 2 )} then f(h) = 
-h/AH 2 . 

For short times, formulae © and 10} are still correct in the presence of 
the force f(h). Indeed, the displacement of the random walker is the sum of a 
random part, of the order of \J~Dt, and a drift part Df(h)t ~ —Dht/AH 2 . The 
former dominates the latter for short times. 



2 Very slow nuclear spin relaxation 

The resonance condition \H z (t) — H x \ < e is satisfied during a time of the order 
of e 2 / D. The flipping probability of a molecular spin during this time is Ae 2 / D. 
If this quantity is of order unity or larger, the nuclear relaxation will be said to 
be 'very slow'. In that case, the proportion n(t) of relaxed spins at time t is 

/OC /" oo 

dhg(h) / dh'p x {hi,h,h\t) (7) 
-oo J —oo 

where p(h\; h, h' , t) is the probability that the random walker, who started from 
h at time 0, is at h' at time t and has reached h% = Hi — H ext at some time t' 
between and t. It is related to p(h, h', t) by 

( p 1 (hi;h,h',t) = p(h,h',t) if (h - h x ){ti - hi) < (a) 

{ (8) 

[ Pi(hi;h,ti ,t) = p(2hi - h,h',t) if (h - hi)[h! - h x ) > (b) 

Relation JS] a) is obvious and JS] b) can easily be derived, or found in 
textbooks UJ. 
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Figure 1: The thick curve shows the hyperfinc field distribution g(AH z ). The 
thin curve shows a possible evolution with time t of the hyperfine field AH Z on 
a particular molecular spin S. 
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The integrand in 10 is very small unless h — hi and h' — hi are of the order 
of \f~Dt or smaller. In that domain the integrand is of the order of l/\/Dt. It 
follows that 5M(t) ~ \f~Dt in agreement with 0. The detailed calculation, 
based on 0, © and 0, yields 

/>oo />oo 

A = 8g(hi)^Dj^ dx dycxp[-(x + y) 2 } (9) 
Jo Jo 

There is no relaxation if \hi\ > AH (if resonances other than at Hi are 
excluded). The maximum value of A (obtained for hi = 0) is of the order of 
1 / y/rjj where 

t h = AH 2 /D (10) 

is the longitudinal nuclear spin relaxation time. The notation Ti has been 
avoided because it usually designates the spin-lattice relaxation time, related 
to spin-phonon interactions, which is extremely long at low temperature. The 
notation has also been avoided because it usually designates the transverse 
nuclear spin relaxation time. 

3 Moderately slow nuclear spin relaxation 

The nuclear spin relaxation will be called moderately slow when the relaxation 
of a molecular, electronic spin requires that (Hi — H z ) vanishes several times, 
but is already almost complete (at resonance, i.e. hi = 0) at t = tjj. 

A qualitative description will be given in the case hi = Hi — H ext = 0. 
As seen from 0, the random walker explores in the time t a held interval 
8h{t) ~ V Dt. The proportion of field values which have been explored at time 
t is thus 

p(t) ~6h(t)/AH . (11) 

This quantity pit) is also the proportion of molecular spins which have a 
chance to relax during the time t. For such a spin, tunneling is allowed during 
a time 

t eff ~et/6h(t) ~ ey/t/D . (12) 

For a spin whose local field lies in the explored region Shit), the relative 
magnetization change, for short times, is At e //, and this evaluation is correct if 
At e // < 1- The total relative magnetization change is for such times 

5M z it)/SM z (oc) =p(t)Xt eff 

or, according to 111(1 and 1)12(1. 

SM z it)/5M z (oo) ~ Xet/AH . (13) 

This formula is different from that written by Miyashita & Saito, who assume 
h = hi at t = 0, so that p(t) is replaced by 1 in IpHfl) . Therefore, they predict 
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5M Z (t) ~ y/t at short time instead of 5M Z (t) ~ i as obtained in ( 1131) . The 
quantity they calculate is actually Xt e ff — XeyVt/D, which is given by their 
formula (3.9) as aJi u>j^y/t/D where a is a constant. Since e ~ Klot / '(g/iBs), 
identification of the two results yields 

A = ojt t x (14) 

where tx is a constant time (related to a). Thus, Miyashita & Saito find that A 
is independent of D and AH. This surprising result is based on the assumption 
of a fairly fast nuclear spin relaxation, and on the statement that "the velocity 
of the field is proportional to \fD n . This point will not be discussed here. 

When t becomes so large that Xt e ft > 1, the depolarization for an 'ex- 
plored' field value is almost total and the average relative demagnetisation is 
SM z (t) /SM z (po) = pit) ~ \fDtj AH — \IWh as in the case of very slow nuclear 
relaxation. This result holds if \hi\ is appreciably smaller than AH. When \hi\ 
approaches AH, the rate of change of the magnetisation can easily be shown to 
decrease to 0. 

4 Long times 

For times t 3> th, the random walker has lost the memory of its initial position 
h. Thus, to quote Miyashita & Saito, "the number of crossings" [of H z (t) 
with Hi per unit time] " is ... constant in time, which causes a constant rate 
relaxation, i.e., the exponential relaxation." Two remarks should be added. 
First, the electronic spin relaxation time r obviously depends on h\ = H\ — H ext . 
Miyashita & Saito made detailed numerical studies, but only in the case of an 
external field tuned at resonance, h\ — 0. When h\ ^ becomes larger than 
the width AH, r obviously goes to oo. The second remark is that, in the case of 
very slow or moderately slow nuclear spin resonance, which is addressed in the 
present note, the electronic spin relaxation is already very strong at t = th , i.e. 
5M (th) /SM (po) is of order unity. This implies that the exponentially relaxing 
part is small. 

5 The case Hi = 

In most of the experiments done so far, Hi = 0, i.e. tunnelling takes place 
between the lowest state of each well. This case is peculiar. 1) The main feature 
is that (if phonons are completely ignored) the system of electronic and nuclear 
spins ignore the temperature T of the crystal, and relax to an equilibrium state 
which generally corresponds to another temperature T e //. Indeed, only a part of 
the Zeeman energy of the molecular spins can be transferred to nuclear degrees 
of freedom. Only for H ext — 0, when the Zeeman energy vanishes, T e // = T. 2) 
If nuclear spin relaxation is not very slow, the above treatment is still acceptable 
for short times, because spin reversal and relaxation are nearly the same thing. 
3) Let the case of very slow nuclear spin relaxation be discussed. Then, when the 
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hyperfine field AH Z crosses the value —H ext , S z follows the field adiabatically 
and changes sign. Thus, if S z (0) — s then S z (t) = s if (h—h\) and (h' — hx) have 
the same sign, while S z (t) = — s if they have a different sign. Thus, the spin 
never loses the memory of its initial state, which is quite unusual in a relaxation 
processs. For short times, formula Q can then be derived from an argument 
similar, but not identical to the above one, and A turns out to be given by 10 
again. For long times, the hyperfine field distribution is probably affected by 
relaxation and depends on time. 

6 Validity of the very slow nuclear relaxation 
scheme 

For a spin able to tunnel between two localized states |— ) and |+), the wave 
function x(t)\+) + y(t)\ — ) satisfies 

i( t ) = ±-x{t)E(+\t)-iu T y(t) ; y(t) = -iu T x{t) + ^y(t)E^ (t) (15) 
in in 

where the unperturbed energies #'(t) satisfy E^(t)-E^(t) = 2gfi B [H z (t)- 
Hi]s. 

Equation simplifies if one introduces the notations u(t) = (1/h) f f dt'E^(t') , 
w(t) = i / t * dt'E^-\t') , x(t) = exp[-iu(t)]X(t) , and y(t) = exp[-iw(t)]Y(t). 
Moreover, the initial condition X(0) — 1 will be assumed, and the time will 
be assumed so short that X(t) may be approximated by X(0) = 1. Then IJ15JI 
yields Y(t) = —iujTexp[—iU(t)] where 

U(t) = {2 9f i B s/h) f dt'[H z {t') - Hi] (16) 

Let the initial value of the local field be tuned so as to allow tunnelling. 
If the tunnelling window remains open until a time of the order of the tunnel 
period I/lot, nuclear relaxation is very slow in the sense defined above. The 
condition for very slow relaxation is thus C/(1/wt) < 1- According to ifTHjl . U(t) 
can be roughly evaluated as the product of (2gfiBs/h) by \J~Dt = AHy/tprn 
and a factor t because of the integration. The condition for very slow relaxation 
reads 

2 gtiB sAH_ < l 

The hyperfine width AH is never smaller than 0.001 Tesla. Thus, i|17fl 
requires large values of th and lot- For instance, a nuclear relaxation time 
th=1 s implies lot > 10 6 s _1 . 

An additional condition is ujtTh > 1, but condition 1)17(1 is probably stronger. 

Usual nuclear spin-lattice relaxation is expected to be extremely slow at 
low temperature. However, measurements by Morello et al. |B] in Mn^ac 
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reveal that th does not increase beyond about 0.01 second. The zero field 
tunnel splitting of the ground doublet does not seem to satisfy the requirement 
lott~h > 1 (Barbara, private communication). Thus, nuclear spin relaxation is 
not very slow and therefore the hyperfine field is not expected to give rise to \ft 
relaxation in zero external field in Mn^ac, but rather to exponential relaxation. 

It is of interest to recall that square root relaxation has been observed in 
Feg on times of the order of a minute. The nuclear spin lattice relaxation 
times T\ which have been reported (for instance when observing negative spin 
temperatures^^) are of the same order of magnitude, but it is not clear whether 
so long nuclear longitudinal relaxation times can be reached in molecular nano- 
magnets. 

7 Real systems: effect of dipole interactions 

In real systems, dipole interactions between molecular spins are present and 
contribute to the local field by an amount which will be called 'dipole field'. 
This contribution adds to the hyperfine field addressed above. Total magnetic 
relaxation requires that the dipole field, as well as the hyperfine field, explores 
the whole allowed region. In Fes, the dipole field created by electronic spins is 
about 10 times as large as the hyperfine field, and therefore 90% of the relaxation 
cannot be explained by the Miyashita-Saito mechanism described above, but 
only by the Prokofiev-Stamp mechanism. 

8 Conclusion 

In the present work it is shown that the \ft behaviour can in principle arise from 
hyperfine interactions alone. This can happen for any value of the external field 
and for any initial state provided it is not the equilibrium state. However, very 
slow nuclear spin relaxation is necessary. In the case of moderately slow nuclear 
spin relaxation, the decay is found to be linear at short times, in contrast with 
the statements of Miyashita & Saito. This discrepancy occurs because Miyashita 
& Saito assumed the initial condition h — hi which is not fulfilled by nuclear 
spins. After some time the decay crosses over to the \fi behaviour and finally 
becomes exponential at long times. The constant A of formula (JTJ, and the 
relaxation rate 1/r which describes the long time behaviour, go to when the 
tuning parameter (H ext — Hi ) becomes large with respect to the hyperfine width 
AH. 

When Hi = and H ext ~ 0, spin relaxation has remarkable features when 
the spin-lattice relaxation time is very long. The spin temperature can then be 
quite different from the lattice temperature. 

The present work contains several shortcomings. The possibility of simul- 
taneous reversal of the electronic spin and a few nuclear spins has been dis- 
regarded. Such processes have been treated by Prokofiev & Stamp^T] in the 
case of superparamagnetic grains. They might be less crucial for molecular 
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nanomagnets since each molecular spin interacts in a significant way with a re- 
stricted number of nuclear spins. On the other hand, the possible dependence 
of D and A with respect to the local field has also been ignored (as in the work 
of Miyashita & Saito). Such a dependence would modify the value I© of the 
coefficient A. 

The present theory may be relevant, for instance, in two cases, i) Mn^ac if 
a transverse field is applied in order to increase lot and to fulfill the condition 
<(T7I) of very slow nuclear spin relaxation, ii) Diluted samples of Fes, where the 
dipole interaction between molecular spins is smaller than the hyperfine field. 
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Wolfgang Wernsdorfer 
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